HUPD-9912 



1 



Neutrino Mixing in the Seesaw Model 

Junya Hashida Takuya MOROZUMi S and Agus Purwanto 



Department of Physics, Hiroshima University, Higashi- Hiroshima 739-8526, Japan 

We study the neutrino mixing matrix (the MNS matrix) in the seesaw modeL By as- 
suming a large mass hierarchy for the heavy right-handed Majorana mass, we show that, in 
the diagonal Majorana base, the MNS matrix is determined by a unitary matrix, S, which 
transforms the neutrino Yukawa matrix, y^, into a triangular form, j/a- The mixing matrix 
of light leptons is VkmS'* , where Vkm = VlJ^Vlv and S'* = Vlv'^ S* . Large mixing may 
occur without fine tuning of the matrix elements of even if the usual KM-like matrix 
Vkm is given by Vkm = 1- This large mixing naturally may satisfy the experimental lower 
bound of the mixing implied by the atmospheric neutrino oscillation. 



§1. Introduction 

Since the observation of neutrino oscillations, 0)'li''& the mixing matrix of the 
lepton sector (the MNS matrixa)) has been discussed in many works. Among the var- 
ious proposed explanations of the light neutrino mass, the seesaw mechanismB-* is an 
attractive scenario. In this paper, we study the MNS matrix in the framework of the 
seesaw model. Specifically, we assume that the gauge group is SU{2) x U(l). Since 
there are three light neutrinos, each one may have its partner, a gauge singlet neu- 
trino with large Majorana mass. Mor eover , neutrino oscillation experiments 



suggest that there is mass hierarchy. r***^l In the context of the seesaw model, the 



mass hierarchy may originate from the Yukawa term and/or Majorana mass term. 
To account for the hierarchy and the mixing, there are two extreme explanations. 
The first explanation is that the hierarchy comes from the Majorana mass term. The 
other explanation is that the heavy right-handed neutrinos have a degenerate mass 
and the hierarchy comes from the Yukawa term. In this paper, we are interested 
in the former case. If the mass hierarchy of neutrinos comes from the hierarchy of 
the Majorana mass term {thni S> mN2 ^ ruNs), the diagonalization of the seesaw 
matrix is carried out in a very different way than that for charged fermions. This 
is what we discuss in this paper. The seesaw mass matrix is a 2Nf x 2Nf matrix 
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(Nf = 3 in our case) and has the simple structure 



M.= . (1) 

Its Nf X Nf submatrix corresponding to the Majorana mass term of light neutrinos is 
zero, and the submatrix corresponding to the Majorana mass term of heavy neutrinos 
can be a real diagonal matrix. The origin of the flavor mixing is the Dirac-type 
Yukawa term denoted by Vv^- Because the seesaw matrix is a symmetric matrix, 

it can be diagonalized by a unitary matrix V as y^Mi.y.ll' Despite the simple 
structure, it is not possible to diagonalize the matrix analytically, because we have 
to treat a 6 x 6 matrix. We perform the approximate diagonalization of the mass 
matrix and obtain the parametrization of V . This produces the mixing matrix of 
neutrino sector. Combining it with the mixing matrix of the charged lepton sector, 
we may determine the flavor mixing in the lepton sector, namely, the MNS matrix. 
The method we employ is very similar to the diagonalization of the seesaw matrix for 
the quark mass in the context of a left-right model. IIP In the model, isosinglet quarks 
with a large mass hierarchy play a role similar to the right-handed heavy neutrino 
in the present case. An approximate diagonalization procedure of the seesaw type 
matrix is developed in Ref. 8). In this paper, we extend that method to the seesaw 
model for the neutrino mass with the same gauge group and higgs as in the standard 
model. 

The paper is organized as follows. In section 2, by demonstrating the procedure 
of the diagonalization in the basis in which the Majorana mass term is diagonalized, 
we explain how the MNS matrix comes out by introducing a triangular matrix. In 
section 3, we study the special case that the Yukawa matrices for charged leptons and 
for neutrinos are simultaneously diagonalized through a biunitary transformation. 
If this were the case for the two Yukawa matrices of up- and down-type quarks, 
the Kobayashi-Maskawa (KM) matrix would be trivial. However, in the seesaw 
model for the neutrino mass, large mixing may still occur. We also discuss the 
phenomenological implications of our analysis. In section 4, we study the large 
mixing situation in the weak basis in which the Yukawa matrices of charged leptons 
and neutrinos are diagonal. In this basis, we study the texture of the Majorana 
mass matrix, which leads to the large mixing, and by doing so, we can understand 
the origin of the large mixing from a different angle. Finally, our conclusions are 
presented in section 5. 

§2. The MNS matrix in the Seesaw model 

We start with the mass terms for lepton sector. Without loss of generality, the 
Majorana mass matrix m^r is assumed to be a real diagonal matrix, rriNij — ^Ni^ij- 
& Writing AI^ in Eq. (1) in the diagonal Majorana mass basis, the mass terms are 

~ ^lepton — ^^Ve'^^^R 
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T V 



rriNi 



•mN2 



(2) 



We focus on the case that the Majorana mass term is much larger than Dirac mass 
term: rriM 3> Uu-^- We are interested in the case that the diagonal Majorana masses 
satisfy the relations rriNi ^ uiNj {i < j), and the rank of the Yukawa matrices ye 
and Ui, is Nj. 

The approximate diagonalization procedure of the seesaw-type matrix goes as 
follows. Mjy is diagonalized by a 2Nf x 2A^£ unitary matrix as V'^M^V = -/Vf^^diagonai- 
Let us parametrize this unitary matrix asEl^ 



V 



S K 
T A 



(3) 



where 5, T and A are Nf x Nf submatrices. The unitarity relation is written in 
terms of the submatrices: 



SS^ + KK^ = 1, 

rrt + AA^ = 1, 

ST^ + KA^ = 0. 



(4) 
(5) 
(6) 



Then M^diagonai is given by 



M.diagonal = ^ M,V = 

( TV^S + S^yuf^T + T^m^T T^yJf^K + S^y^f^A + T^m^A 
\ (tV^.K + S^y^^^A + T^mr^Af A^yJ^K + KTy,^A + A^m^A 



'V2 



V2 



We can easily see that the submatrices 



T^yJ^S + S' y^^T + T' niNT. 



V2 



V2 



A^yJ-^K + K'^y^-^A + A^tunA 
must be diagonal, and the off-diagonal submatrix becomes zero: 



(7) 

(8) 
(9) 



T^y/^K + S'^yu^A + T^tunA = 0. 



(10) 
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In order for the O(mjv) part of Eq. (9) to be diagonal, we must have A = 1 , Ell 
because m^r is diagonal and its diagonal elements are not degenerate. Then Eq. (5) 
implies that the leading order of T is 0(;^^). Equation (6) means K = 0{:^)S. 
Therefore, to leading order, 5*^5 = 1. In order for the 0{v) contribution of the RHS 
of Eq. (10) to be canceled, we can neglect the contribution of the first term. Then 
the relation T = — y^S is obtained. K is also determined through Eq. (6). 
The above consideration leads to the following submatrices: 



A = l-0 



niN 



SS^ = 1-0 



rriN 



y^S, K = 



(11) 

Below, we neglect the deviation from unity of A and SS'^ . Now, the seesaw mass 
matrix is written as 



-QT To 







mN + 0{v) 

At this stage, the unitary matrix is parametrized as 



V 



1 



\/2»njv 



T 




(12) 



(13) 



Further, 5 must be chosen so that the —S^^yu-^^Vu S is diagonal. For any matrix 
y of rank Nf, we can find a unitary matrix which transforms y into a triangular 
form.H^ We choose such matrix as S'^: 



S'^yv 




= VA- 



(14) 



Here the diagonal elements are real, and the off-diagonal elements are complex in 
general. In this form, it is noted that the lightest neutrino Ve only interacts with 
A''!, couples to both Ni and N2, and Ur couples to all heavy neutrinos Ni, N2 
and A3, as shown in Fig. 1. By choosing 5 in this way, the seesaw mass matrix is 
approximately diagonalized. This is shown by writing submatrix —S'^yyi^yJ^S: 



rriNi 

yiy2i 

mpfi 

yiyai 

mjvi 



y2 



"IJV2 

y2iy3i 

mjvi 



yiy2i 

+ 



yiy3i 

mjvi 



y2i 



*' We can also choose another A, 



A 



mjvi 
_|_ i/2i/32 



1 





y2iy3i 



+ 



y2y32 

mN2 



y3 



+ 



y32 

m,N2 



+ 



y3i 

mjvi 



.(15) 



This choice changes only K as 



K = ~ST^ A, and leaves the Majorana mass term in a diagonal form, but with a reversed sequence of 
mass elements, i.e. diag(mjvi, mjv2, 'tijvs) — > diag(mjv3, mjv2, ?tijvi). In our approximation neither 
S nor T changes. This fact in turn implies a light neutrino submatrix [Eq. (8)] of the same form as 
in the diagonal A case. 
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Because mjvi 3> m^2 ^ 't^n^^ the off-diagonal elements are much smaller than 
the difference of the diagonal elements. Therefore the mass matrix is regarded as 
approximately diagonal: 



/ -ml 



\ 



2 2 



2mjv3 



mm 



mN2 



V 



(16) 



mm j 



The above matrices show that the mass of is determined by the mass of the 
heaviest singlet neutrino, A'^i, and is not affected by the presence of the other lighter 
singlet neutrinos, because v^. couples to only A''!. Moreover, though i/^ couples to 
both A'^i and N2, its mass is mainly determined by the mass of A''2, because the mass 
of A*"! is much larger than that of A^2- Similarly, the mass of the Nj, mainly depends 
on the mass of A^a, due to the relation mjsfi ^ 171^2 ^ n^Ns. If such a mass hierarchy 
of the three heavy neutrinos is admitted, this procedure is a powerful method for 
diagonalizing the seesaw-type mass matrix. 

Then, the unitary matrix V can be rewritten as 



V 




1 



-J/A 



^ 1 



(17) 



The MNS matrix is a 3 x 6 submatrix of V multiplied by a 3 X 3 unitary matrix for 
diagonalization of the Yukawa term of the charged Icpton, 



VMNS = iVLjS*){l,yA 



(18) 



Here Vie comes from the diagonalization of the Yukawa term of the charged lepton 
and is defined by the equations 



e\ = VReCR, 
VeD = VLjyeVRe, 



(19) 
(20) 
(21) 



where peD is a real diagonal matrix. We introduce a unitary matrix which transforms 
the Yukawa term of neutrinos to the diagonal form yuD- 

yu = VluVudVrJ 

= S*yA. (22) 

The second line of Eq. (22) comes from Eq. (14). Therefore we can relate t/a and 
yi^D through the equation 



yuD = VLjS*yAVR„ 
= S'*yAVR., 



(23) 



5 



where we introduce another unitary matrix S' , 



S'* = VlJS*. 



(24) 



S' represents the difference between the unitary matrix which transforms yi, into 
the triangular form and the unitary matrix for the diagonaUzation by a biunitary 
transformation. 

Finally, we can write the leptonic charged current as 



Jlh = {eLlJ-LTL)jnVMNS 



( VLX 

1 

Nr2 



\ 



(25) 



with 



VMNS = {VLjS*)il,yA 



= iVKMS'*)il,yA 



V 



\/2mjq' 



and 



Vkm = VlJVi 



Lv 



(26) 



(27) 



Here vl and Njf are mass eigenstates of neutrinos. Prom the first definition of Vmns 
from Eq. (25), we can see that the mixing of the neutrino sector is determined by 
the unitary matrix which transforms the neutrino Yukawa term into the triangular 
form. In the second definition of Vmns from Eq. (25), we introduce the usual KM- 
like matrix and rewrite the matrix Vlc^ S* as {VLe'^VLv){VLj S*) = VkmS'* . By 
rewriting the matrix in this way, we can divide the Icpton mixing angles into two 
parts, one which comes from the usual KM-like part, denoted by Vkm, and another 
(S") which comes from the difference between the unitary matrix Vlv for biunitary 
transformation and the unitary matrix S for transformation into the triangular form. 
The advantage of the latter decomposition is that both Vkm and S' are independent 
of the choice of the weak basis. 

We now summarize the results obtained in this section. The MNS matrix is a 
product of the mixing matrix of the charged lepton sector and that of the neutrino 
sector. We find that the unitary matrix denoted by Vlc^ S* determines the mixing 
matrix among light leptons. Here S denotes the unitary matrix which transforms 
the neutrino Yukawa term into a triangular matrix. Vlc is the unitary matrix for 
the diagonalization of the charged Icpton Yukawa term. This is proved in the weak 
basis in which Majorana mass matrix of heavy neutrinos is real diagonal and also 
under the assumption that the heavy Majorana neutrinos have greatly differing mass 
scales. 
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§3. Large mixing 



The MNS matrix obtained in the previous section is divided into two parts 
unitary matrices S' and Vkm- We note that S' is determined by yu and niN in the 
diagonal Majorana basis, while Vkm is determined by yuVu^ and VeVe'^ ■ Therefore S' 
is determined by the mass matrix of the neutrino sector, while Vkm is determined 
by the Yukawa terms only. In order for Vk m to be a non-trivial matrix, the following 
condition must be satisfied: 

Vkm^i ^ [yeyJ,y^yJ]^o. (28) 

On the other hand, the condition for S' being non-trivial is 

S'^l ^ {yAyA^h^O for (i^j), (29) 

in the diagonal Majorana basis. Equation (29) comes from Eq. (23), because 

S'*yAyA^S''^ = yuD^- (30) 

The mixing angle is the product of S' and Vkm- We note that for a given t/a we 
may determine S' . Equation (29) and ( |30| ) imply that the off-diagonal elements in 
the triangular matrix induce the mixing in S' . If the off-diagonal elements of 
vanish, S' = 1. However, since Vkm is arbitrary, in general, lepton mixing angles 
are also arbitrary. Even if there is large mixing in 5", this can be canceled by Vkm- 
There may also be the case that S' is the unit matrix and there is large mixing 
in Vkm- Here we focus on the case that Vkm = 1- Equivalently, we impose the 
commutation relation 

[yeye\yuyu] = ^- (31) 

If this is the case, S' determines the mixing angles. 

Now we show that large mixing occurs without fine tuning of the matrix elements 
of y/\. We illustrate this by considering the two flavor case {Nj = 2) and apply Eq. 
(31) to the explanation of the large mixing of atmospheric neutrinos. In the two flavor 
case, yAl/A' ^' ^"^^ approximately diagonalized light neutrino mass matrix my 
are respectively expressed as 

y^yi = { , 2 \ , (32) 



s' 



cos 9 
sin0e~* 



sin Oe^'t' 
cos 9 



(33) 



with 



tan 26* 



2|y32|y2 



-^2 + 1^32^ + 2/3' 



arg(y; 



'32 j 



(34) 
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and 



y2y32 \ ^2 



I y2y32 ys , yz2 I o ' v / 

mjV2 mjv3 niJV2 

where m7V2 ^ ^7V3 is assumed. In order that Eq. (35) can be regarded as an ap- 
proximately diagonal matrix, the sufficient conditions on the elements of the Yukawa 
term are 

0{v2) = 0(y3) ^ y, (36) 

O {^^^ < 1- (37) 

The first condition is that the hierarchy of the light neutrino masses comes from the 
hierarchy of the Majorana mass term, rn-7V2 ^ t^NZ- The second condition suppresses 
further mixing from the form of Eq. (35). 

In Fig. 2, we show sin^ 26 as a function of ■'^^ by changing the ratio ^ = 
For r = 1, maximal mixing occurs at infinitesimally small |y32| = e, and the mixing 
decreases as |y32| increases. For r < 1, maximal mixing occurs at non-zero 7/32 
because the diagonal elements are degenerate for |y32| 7^ 0. For r > 1, the mixing 
angle is suppressed because the diagonal elements cannot be degenerate. For r / 1, 
small mixing is also possible for |y32| <C {ys — y2\- In Fig. 3, we constrain the ratios 
of the Yukawa couplings by using the result obtained by the Super-Kamiokande, Eq. 
(41). The allowed range is between two curves in the parameter space, (r, ■^^^)- We 
can see that the experimental constraint can be easily satisfied for the range of the 
parameters in Eqs. (36) and (37). The curves are obtained with the equation 

(S)^(^~t^) =^ + to?2^' ^^^^ 

where sin^ 29 takes the minimum and the maximum values of the experimental con- 
straint, Eq. (41). The qualitative features of Figs. 2 and 3 can be understood by 
using the approximate formulae for tan 29 under the conditions set by Eq. (36) and 
Eq. (37). From Eq. (34) we have 

r 3^>2 for |y32p>yi-yi, 
tan20= 'I ^' (39) 

I S |y32p <yi-2/i. 

The first line corresponds to the case that the diagonal elements are very degenerate 
and the second line to the case that they are not degenerate. The first case of Eq. 
(39) leads to 

sin^ 29 > 0.8. (40) 

This can be compared with the result for atmospheric neutrinos. The Super-Kamiokande 
collaboration & has reported the range 

sin^ 29^tm > 0.8. (41) 
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Moreover, the second line in Eq. (39) means that large mixing occurs for |y32| ~ 

lys - 2/21- 

Finally, we comment on the mass scale rriNs of the lightest gauge singlet neu- 
trino. According to the Super-Kamiokande result, the mass-squared difference of 
the pL and r neutrinos is of order 10~^[eV^]. Assuming that this corresponds to the 
mass squared of the r neutrino and that the Yukawa coupling is the same order of 
magnitude as the gauge coupling, we obtain 

m^r - 42 
mm 

Therefore, the mass of the hghtest gauge singlet neutrino may be 0(10^^) [GeV]. 



§4. Texture of the Major ana mass 

In the previous section, we obtained large mixing under the additional constraint 
that the Yukawa matrices of charged leptons, yevl, and neutrinos, Uuvt^ ^'I'c simul- 
taneously diagonalized, i.e. Vl^ = Vl^ [see Eq. (31)]. The origin of the mixing may 
be traced to the presence of the off-diagonal element of the triangular matrix in the 
diagonal Majorana mass basis. It is useful to study the same situation from a dif- 
ferent viewpoint. In this section, we study the large mixing case in a different basis. 
We can transform from the diagonal Majorana mass basis to the diagonal basis of 
Yukawa terms for charged leptons and neutrinos. If we further impose the condition 
Eq. (31), the transformation becomes a transformation of the weak basis. Therefore 
Vmns becomes the 3x3 unit matrix. In the new basis, the Majorana mass term 
is not diagonal, and we may investigate the texture of the Majorana matrix that 
induces large mixing. 

In this basis, the mass term (2) is rewritten as 



V 



y>^D^o 




where 



= VlJuI (44) 
K = VrJnI (45) 



and the leptonic charged current is 



'L3 

Up to this point, we have not imposed the condition Eq. (31). If this condition is 
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imposed, namely Vkm = 1; the charged current becomes 



Jlh = (eLjlLfLhn I ^^'12 1 • (47) 

Vru'^^tt-nVru (= "?,jv) is the Majorana mass matrix in the basis of the diagonal 
Yukawa term. Because of Eq. (23), Vjii, is determined by ua^Va as 

Vud'^ = VR^^VA^yAVRu. (48) 

Now we give the form of the Majorana mass matrix that induces large mixing 
in the two flavor case {Nf = 2). Because ua^Va is expressed as 

V 2/322/3 2/3 / 



we can parametrize Vr^, as 



cos a — sm a e^^ 



= I ^, I ' (50) 



sma e cos a 



with 



tan 2a = '^[^"^lo'' 2' = -arg(y32)- (51) 

+ 1^32 r -y3^ 



Therefore, the Majorana mass matrix is written as 

_( mAr2 cos^ a + miV3 sin^ a e~^*'^ — |(mjv2e*'^ — JTijvse"*'^) sin 2a 
i — ^("7,jv2e*^ — 'mNse~'^^) sin2a mAr2 sin^ a e^*'' + mAr3 cos^ a 



(52) 



In Fig. 4, we display the sin^ 2a as a function -l^J- by changing the ratio f = 
By comparing Figs. 3 and 4, we can see that maximal mixing of the mixing angle 
9 and maximal mixing of a take place at the same time only when the diagonal 
elements of yA are exactly degenerate (r = 1). In this case, by setting sin^ 2a ~ 1, 
the Majorana mass matrix can be expressed as 



mm 1 + 66-2^/3 _(i _ ee-2^/5)e^/5 . 



where e = mm/mm (^ !)• Prom Eq. (53), the texture which leads to large mixing 
has the characteristic feature that the off-diagonal and diagonal elements are the 
same order of magnitude and the Majorana mass matrix is close to a democratic- 
type matrix. There are small corrections to the democratic- type matrix, which lead 
to the non-zero mass of the heavy neutrino -/V3. 
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§5. Conclusions and discussion 

We have studied the MNS matrix in the seesaw model. We proposed an approx- 
imate procedure for diagonalization of the seesaw matrix. If we start with a real 
diagonal Major ana mass matrix and assume a hierarchy for the masses of the heavy 
neutrinos, we find that the unitary matrix S that transforms the Yukawa term into 
a triangular form determines the mass eigenstate of light neutrinos up to a small 
mixing of the heavy neutrinos. The triangular form has to be arranged in such a 
way that the lighter SU{2) doublet neutrinos couple to the heavier singlet neutrinos. 
The mixing matrix for light leptons is V^^S*, i.e., the product of the unitary matrix 
Vie required by diagonalization of the Yukawa matrices of charged leptons, yevl, 
and the unitary matrix S* to transform the Yukawa matrix of neutrinos, yu, into 
the triangular form y/\. By inserting the unitary matrix for the diagonalization of 
the Yukawa matrices of the neutrinos yuyt, the mixing matrix can be rewritten as 
VkmS'* ■ Vkm is a usual KM-like matrix and denotes the difference between the uni- 
tary matrices for the diagonalization of Yukawa terms, and S' denotes the difference 
between the unitary matrix for the transformation into the triangular form and the 
unitary matrix for the diagonalization of y^- Therefore S' can be determined by the 
structure of the mass matrix of the neutrino sector only. We find that for a given y^x , 
we can determine S' . We also find that the presence of the off-diagonal element of t/a 
induces the mixing in S'. A large mixing angle is possible and it was explicitly given 
for the two- flavor case. Since Vkm is arbitrary, we impose the additional condition 
Eq. (31) so that S' determines the mixing matrix. By doing so, we compare the 
obtained mixing angle with the atmospheric neutrino anomaly, which was reported 
by the Super-Kamiokande collaboration. Because the mixing angle is a function of 
the elements of y/\ , we can find the region of the elements which is compatible with 
the Super-Kamiokande results. We would like to comment that the origin of the 
commutation relation of the Yukawa terms may be the unification of the quark and 
lepton. Because the KM matrix is close to the unit matrix, this may suggest that 
up and down quark Yukawa terms may approximately satisfy similar commutation 
relations. Such a relation may also hold for the lepton sector in unified theories of 
quarks and leptons. The phenomenological study of such systems can be extended 
to the three flavor case and may be applied to the solar neutrino problem. We also 
note that Smirnov0) studied mixing in the seesaw model. Because our approach 
differs from his, and the transformation to the triangular form is essential in our 
approach, the relation between our work and his must be carefully examined. We 
will study these issues in future publications. 
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Fig. 1. In the form of j/a, only couples to A^i, while couples to both A''i and A^2, and Vt couples 
to all heavy neutrinos, A''i, A^2 and N3. The mass of i/e depends only on mjvi. Moreover, due 
to the mass hierarchy mjvi ^ mN2 3> mjva, and j/t are determined by the masses of mN2 
and niNS, respectively. 




Fie. 2. sin^ 20 as a function of for various r = — . The three curves correspond to r = 1 

^ y2 V2 ^ 

(solid line), r = 1.2 (dotted line), and r — 0.8 (dot-dashed line). The region between the two 
solid straight lines is allowed by the Super-Kamiokande experiment on atmospheric neutrinos, 
0.8 < sin^ 2eatm < 1. 
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Fig. 3. The region in the parameter space -^^^) between the dashed lines is allowed from the 
Super-Kamiokande experiment on atmospheric neutrinos, sin^ 26atm > 0.8. 
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Fig. 4. sin 2a as a function of ■'^^^ for various ^ ~ The three curves correspond to r = 1 (solid 
line), r — 1.2 (dotted line), and r = 0.8 (dot-dashed line). 
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